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' Abstract 

^SJ ■ Using the Cornwall-Jackiw-Tomboulis effective potential for composite 

> , 

Q^ , operators we compute the QCD vacuum energy as a function of the dynam- 

a^ 

^S| ' ical quark and gluon propagators, which are related to their respective con- 

(^ , densates as predicted by the operator product expansion. The identification 

o\ 

Q^ ■ of this result to the vacuum energy obtained from the trace of the energy- 



Ph! 
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momentum tensor allow us to study the gluon self-energy, verifying that it 



Q_i' is fairly represented in the ultraviolet by the asymptotic behavior predicted 

,i^ . by the operator product expansion, and in the infrared it is frozen at its 



asymptotic value at one scale of the order of the dynamical gluon mass. We 
also discuss the implications of this identity for heavy and light quarks. For 
heavy quarks we recover, through the vacuum energy calculation, the relation 
nif < Tpftpf >~ — Y2 i^G^'^G^u) obtained many years ago with QCD sum 
rules. 
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1 Introduction 

The QCD vacuum energy can be computed through the use of the effective potential 
for composite operators [|l| . This vacuum energy is a function of the masses of the 
theory |]T], |^. For quarks the effective mass is a sum of the running mass and the 
dynamical one and the quark condensate is an order parameter for chiral symmetry 
breaking (see, e.g., 0]). According to the operator product expansion (OPE) [|], Q 
the dynamical mass is a function of the quark condensate. Therefore, the vacuum 
energy involves the bare quark masses and the quark condensate. 

The concept of a dynamical gluon mass is crucial to calculate the QCD 
vacuum energy. (We would like to emphasize that the presence of a dynamically 
generated mass does not mean that gluons can be considered as massive asymptotic 
states like dynamically generated quark mass does not mean that quarks can be 
observed as massive asymptotic states. Why quarks and gluons are not observed 
as free states is the well known problem of confinement). As discussed at length in 
Ref. 0, the vacuum energy is finite as long as the gluon mass decreases at large 
momentum, and it vanishes when the gluon mass is zero. The infrared divergences 
that plagued the effective potential calculations are also absent for a finite dynamical 
gluon mass. Finally, the gluon mass is related to the gluon condensate, and a precise 
relation, as in the case of the dynamical quark mass, has been obtained through the 
OPE P] . It is also known that the gluon contribution dominates the vacuum energy. 
However, not much is known about the infrared behavior of the dynamical gluon 
mass. 

On the other hand the vacuum energy can also be obtained through the trace of 
the energy momentum tensor and is a function of the quark and gluon condensates 
and the bare quark mass. Equalizing these two expressions we obtain information 
about the dynamical masses of the theory and the condensates. In this work we 
use two different calculations of the vacuum energy density in order to determine 



the behavior of the dynamical masses, particularly of the gluon dynamical mass. 
We show that it is fairly represented in the ultraviolet by the asymptotic behavior 
predicted by the OPE, and in the infrared it is equal to its value at a scale of the 
order of the dynamical gluon mass. Within this approach we are also able to obtain 
a relation between the quark and gluon condensates in the heavy quark limit. It 
is oportune to recall that many years ago Shifman et al. introduced the notion 
of a nonvanishing gluon condensate |^. Making use of the OPE they started the 
successful program of QCD sum rules, where observable quantities can be written 
as a function of quark and gluon condensates, and one of their main results is the 
relation between quark an gluon condensates valid in the heavy quark limit |10[|: 



mf<i;fi,j>r^-^^^G^-^G,^y (1) 

This relation emerges naturally in our calculation of the vacuum energy. We also 
discuss what can be learned from our approach in the chiral and light fermion limits. 
In Section 2 we setup the equations for the vacuum energy. Section 3 contains a 
discussion of our computational ansatze for the dynamical masses, and Section 4 
contains the vacuum energy calculation. In the last section we discuss the method 
and present the conclusions. 

2 The effective potential for composite operators 
and the QCD vacuum energy 

For a non-Abelian gauge theory the effective potential has the form [||] 
V{S,D,G) = -ijAl-Jr{\nSo'S-So'S+l) 

+V2{S,D,G), (2) 



where S, D and G are the complete propagators of fermions, gauge bosons, and 
Faddeev-Popov ghosts, respectively; Sq, Dq, and Gq the corresponding bare prop- 
agators. ¥2(8, D,G) is the sum of all two-particle irreducible vacuum diagrams, 
depicted in Fig.l, and the equations 

Js'Id'Jg'' ^ ^ 

are the Schwinger- Dyson equations (SDE) for fermions, gauge bosons, and ghosts. 
We can represent ^(5", D, G) analytically by 

iV2{S,D,G) = -^Tr{TSTSD)-^Tr{FGFGD) 

+ irr(r(3)Dr(3)DD) + -Tr{T^*^DD) (4) 

where F, F, T^^^ and F*^^) are respectively the proper vertex of fermions, ghosts, 
trilinear and quartic gauge boson couplings 0. In Eq.(|D we have not written the 
gauge and Lorentz indices, as well as the momentum integrals. 

The complete gauge boson propagator D is related to the free propagator by 

D-' = D,'-U^, (5) 

where 11^ is the gluon polarization tensor, which is obtained from Eq.(^ and Eq.(^), 
and described by 

n^ = F^F^ + FGFG - -T'^^'^DT'^^^D - -T^^^D. (6) 

2 2 ^ ^ 

The diagrams contributing to 11"^ are shown in Fig. 2, and this self-energy is the one 
that admits a massive solution as described by Cornwall 0. The Landau gauge 
expression for the complete gluon propagator is 

The complete fermion propagator S is related to the free propagator by 



^''V) = -I^^^(^''^-^)- (7) 



S-' = So' - S, (8) 



where Sq = i/{^ — mo) {rriQ is the bare quark mass) and S is the Schwinger-Dyson 
equation for the quark self-energy, which is given by 

E = -^STD, (9) 

where 7 is the bare vertex. 

The vacuum energy density is given by the effective potential calculated at min- 
imum subtracted by its perturbative part, which does not contribute to dynamical 
mass generation 0, |^, 0] 

(17) = Vmin{S, D, G) — Vmin{Sp, Dp, Gp), (10) 

where Sp is the perturbative counterpart of S, etc... 

Vmin is obtained substituting the solutions of Eq.(^ into Eq.(^. It is oportune to 
recall that we shall compute the vacuum energy by using ansatze for the dynamical 
masses. Although these ansatze are consistent with the OPE they are simple ap- 
proximations to the full solution of the Schwinger-Dyson equations. However, since 
we calculate the potential at its stationary point, the result depends significantly 
less on the approximations used, as was noticed for the vacuum energy calculation in 
the case of fermions W^ . If we use the corresponding Schwinger-Dyson equations. 



Vmin can be divided in two parts (as will become clear in the following) 

V ■ = V^ ■ + V^ ■ fill 

ymin ^ mm ^ ^ mm) \^^) 

where the labels / and g indicate the fermionic and gluonic parts of the vacuum 
energy. The quark contribution to the effective potential obtained by substituting 
the solution of 5V/5S into Eq.(^ and (in Euchdean space with P^ = — p^) is equal 



too 






(27r)^ 



— ln(' 



(12) 



'P'^ + ml' P^ + T."^ 
where N is the number of colors (quarks are in the fundamental representation of 

SU(N)). 



To compute V^j„ we make an approximation which considerably simphfies the 
calculation. This approximation frequently performed to solve the gluon polariza- 
tion tensor is due to Mandelstam [jl2[ and consists in neglecting the ghosts and the 
diagram with quartic coupling. The neglect of ghosts diagrams was shown to be rea- 
sonable, because their contribution is numerically small [Q. These approximations, 
with the use of the Landau gauge, were shown to be satisfactory in the lengthy and 
detailed work of Brown and Pennington [Q. Although it may seems rough this 
approximation has been verified to be phenomenologically consistent in Ref. [0 as 



long as the vacuum polarization tensor involves a mass scale. Furthermore, as can 
be easily verified, by expanding the integrand in Eq. ( ]T5| ) in powers of 11 we obtain 
a result similar to the Cornwall determination of the vacuum energy in the case of 
a pure gluon theory 0]. Therefore, only the diagram with trilinear coupling should 
be considered in Eq.(H), i.e. the gluon polarization tensor is given by 



2 



(13) 



With this approximation the two-loop gluon contribution in Eq(^) is reduced to 
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Finally, we obtain |n 



(2^ 



n , / n 



n^ 



3P2(p2 + n) 



(14) 



(15) 



where all the quantities are in Euclidean space and A^ = 3 for QCD. 
Thus, according to Eq.(lO), the vacuum energy is 



(16) 



where {Vtf) is (see Eqs. (10) and (12)) 
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(17) 



where m'j = m'j{P'^) is the running quark mass, which is the perturbative solu- 
tion of the fermionic Schwinger-Dyson equation with bare quark mass rriQ. Note 
that Eq.(p!7D has been written for a unique fermion of mass rrif. If we consider nf 
fermions of the same mass it should be obviously multiplied by Uf. There is only 
one divergence in (Qf) that is associated to the bare quark mass renormalization 
as will be discussed afterwards. For gluons with dynamically generated mass the 
gluon contribution to the vacuum energy is free of divergences [0, |T^. As long as 
we have appropriate expressions for E and 11 the calculation of (Q) is an easy task. 

3 The self-energy of quarks and gluons 

To compute the vacuum energy we need reasonable ansatze for the quark and 
gluon self-energies. In the case of quarks the high energy behavior of the dynamical 
mass given by the OPE is P, || 

m.yniP') - ^^^^^/(^^)^^(lnPVAy-^-, (18) 

where g'^{P'^) is the running coupling constant 

where b = 2Att'^ / {UN — 2nf). In the framework of SDK the same behavior was 



found in Ref. [0. As well known the asymptotic behavior of the running quark 
mass is the following: 

-r-(P')-'M.'){^) . (20) 

where Aq is the QCD scale, 7™ = 9C2/(11A^ - 2n/), C2 = {N^ - l)/2N is the 
quadratic Casimir of the fundamental representation, and mflfi"^) is the running 
quark mass at the scale fi^ (a mass scale where perturbative QCD can be safely 
applied), and is related to the bare mass by 

m;(/i)=mo(A)Z„i(^,A), (21) 



where A is a ultraviolet cutoff and Z^^{fi,A) is a renormalization constant for the 
bare mass term 

A discussion of the renormalization procedure applied to the effective potential for 
composite operators can be found in Ref. [|16] . 



Eqs.(|T8|) and (^) involve the high momentum part of the self-energy. To estab- 
lish our ansatz for S for all P^, we recall a few points on the infrared behavior of 



the quark mass. There are several arguments in favor of the relation [1^ 

S(p2 ^ 0) ^ {-ip'^y^^ (23) 

for massless or very light quarks. For heavy quarks we have 

S(p2^0)^m^(p2)|^,^^,. (24) 

These equations provide an idea about the ansatz we should use for the quark 
self-energy. Actually, there is a full numerical solution of the quark Schwinger- 



Dyson equations with bare massive quarks performed by Maris and Roberts [|TS 
and they verify that for very heavy quarks the self-energy is given by Eq.(P^ at low 
momentum and behaves as Eq . (|20|) for a momentum scale larger than the current 



mass, i.e. the contribution of dynamically generated mass to the total mass is 
negligible for heavy quarks. For light quarks {mf < (—%jj^) ) the self energy is a 
sum of the dynamical mass which dominates in the low momentum region but the 
asymptotic behavior of self-energy is still given by Eq. ([20|) . Only at the chiral limit 
S is fully described by the dynamical mass function. This behavior does not change 
when we consider gluons with a dynamically generated mass [l^ , although we obtain 
smaller values for the dynamical quark masses (at P^ — > 0) in this case pO| , pl[| . For 
large gluon masses the signal of chiral symmetry breaking may disappear when we 



solve the Schwinger- Dyson equation in the lowest order |21[], but we consider this as 



a failure of the method for solving gap equations ||2^ at the lowest order and still 
assume that Eqs.(|l3) and (|23|) are valid. 



Our ansatz based on the solutions of Ref. ITSll has the form 



j:{P^)=mf{P^)+mdyn{P^), 



(25) 



which is a sum of the bare and dynamical mass. According to the numerical solutions 



of Ref. |T^ and in agreement with the ansatz already used in Ref. [jTT], [T^ 
have the following expressions: a) In the chiral limit (c/) 



Ki{P^) = m^yr^iP 



we can 



(26) 



where 



mdyn{P^) = rf^ 



^z 



Av' - P' 



1 (\nP^/klV"-~^ , „ „■ 



(27) 



p2 vinWA^y 

where ?] is a mass scale which in principle can be related to the quark condenate, 
since in QCD, according to Eq.(|l^), at some large mass scale (yU >> Ag) we must 
have, r]^ oc (3(A^^ — l)g'^{ii^)/8N'^) (—tjjtp) ■ The above ansatz is equal to the 
one used in Ref. |[l^, where the full effective potential was computed and r] was 
considered a variational parameter to be phenomenologicaly adjusted. In our case 
we calculate the minimum of energy. Therefore, our ansatz must approach as close 
as possible the actual solution of the Schwinger- Dyson equation. Consequently, to 
be also compatible with the low energy phenomenology (Eq. (|23D ), which sets a scale 
for all the chiral observables, we set 

'^' (28) 



b) In the case of a light fermion (//) ( i.e. rrif < rj) we define 

/luPVA^^-^- 



mif{P ) = mf{7] ) 



e{r]' - P' 



e{P' - T]' 



\\nriyAlJ 
which has the correct low and high energy behavior, and use the ansatz 



(29) 



^if 



mif{P^) + mdyn{P^)- 



(30) 



The numerical solutions of the Schwinger- Dyson equations |]18| show that when the 
bare mass is much smaller than the dynamical mass, S(-P^) starts falling as 1/P^ 
just after P^ k. rf^ changing to the logarithmic behavior when rf' jP'^ becomes of 
the order of the bare mass. Eq.(pOD reproduces this behavior and can reasonably 
represent the masses of the quarks tt, d, and s with a suitable choice of bare masses, 
c) For heavy fermions (/i/) we introduce 



m,AP') = nif 






,2 o2 



(31) 



where the basic difference with Eq. ( ^91) is the choice of arguments of the Q functions. 
From the results of Ref. [1^ we could be tempted to set S/j/(P^) = m/^/, which is 
of course the leading contribution to S(P^) for very heavy quarks. However, notice 
that Eq . ([T7|) would give (fij) = with this approximation, and we must stick to 
Eq. (|25|) to obtain the correct result as will be discussed in the next section. 

Eqs. (^) -(|3TD give a quite reasonable fit for the full numerical solutions of the 
quark Schwinger-Dyson equations. The Q functions separate the regions where tech- 
niques like the OPE and renormalization group are reliable from the ones where 
nonperturbative effects are present and we can at most rely on phenomenological 
models. Being a variational parameter the scale ri can be adjusted in order to fit 



observable chiral parameters [16|. Moreover, the freezing of the coupling constant at 



the scale r] is also implicit in Eq. (|27|) . Although there are arguments in the literature 



showing that this freezing may already occur at the gluon mass scale P0| , p^ , we 
do not expect substantial differences in our calculation due to this approximation 
as will be discussed in the conclusion. 

We can now turn to the ansatz for the gluon polarization tensor. The behavior 
of n away from the high momentum region is much less known and is more con- 
troversial. Nonperturbative solutions for the gluon propagator have been searched 
for a long time, but only recently it became more common to talk about a possi- 



ble infrared finite behavior of this propagator. As discussed in Ref. |T^ there are 
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several forms proposed in the literature for the gluon polarization tensor in the in- 
frared, which appear due to different approximations performed when solving the 
Schwinger-Dyson equation for the gluon propagator, and there is an indication that 



the vacuum energy is minimized for gluons with a dynamically generated mass [|T^] . 
A strong support for this possibility comes from the gluon propagator simulation 
in lattice QCD p3[. These results may be considered preliminary and large lattice 
studies must be pursued, but they definitively show signals of a gluon mass scale 
up to reasonably low momenta [^. We shall assume here that gluons have indeed 
a dynamical mass as predicted by Cornwall [^ many years ago, and to obtain an 
ansatz for the gluon polarization tensor we are guided by the OPE again. The OPE 
can teach us only about the high energy behavior of the dynamical mass which is 
given by ^ 



^oMn - - 3^;-:i)'- p."' - (32) 

For the infrared we assume that 11 freezes at the value predicted by Eq. ( P^ at some 
scale which should be determined. This freezing is consistent with the prediction of 
Ref. of a constant gluon mass in the infrared. Therefore, we assume 

n(P^) = f^pixf^l - P') + ^,e{p' - x/i?) (33) 

where 

Notice that x i^ Eq.(^) is a variational parameter to be determined in our cal- 
culation. It defines the scale that separates the perturbative and nonperturbative 



regions ||2^. This choice of ansatz is fully based on the result of the operator product 
expansion for the gluon propagator ||^ and on the phenomenological estimates of the 
dynamical gluon mass ^ |T^ . Note that this is an ansatz in the case of a pure gluon 



theory. Of course, fermions modify this result. There are two types of fermionic 
contributions to the dynamical gluon mass. One is proportional to the gluon con- 
densate but appears only at the next order in the coupling constant, i.e. the gluon 

11 



mass is going to be modified by a factor (1 + asUf...). Tlierefore, for lieavy quarks 
their contribution is small due to the smallness of the coupling constant at the scale 
of heavy quarks. However, for light fermions this is no longer true. Unfortunatelly 
there is not enough information on the infrared behavior of the gluon propagator 
to do better than this. Therefore, we are forced to consider phenomenologically the 
contribution of fermions. The other contribution to the gluon mass is proportional 
to the fermion condensate and its mass mj (ipip)^ which could be important for 
heavy quarks but as we checked is numerically small compared to the contribution 
proportional to the gluon condensate and pure gluon loop. Eqs. (]25|) and (|33| ) are 



the basic ingredients of the vacuum energy computation that we present in the next 
section. 

4 The gluon self-energy and a relation between 
the quark and gluon condensates 

As discussed in the beginning our main goals are to study the infrared behavior of 
the gluon propagator and obtain a relation between the quark and gluon condensates 
like Eq. (|l|) through the vacuum energy analysis. For this we recall that the vacuum 
expectation value of the trace of the energy momentum tensor of QCD is [^] 

(©mm) = ^ {G,uGn - E^/ {-^^) (1 + 7rn), (35) 

^9 f 



where the perturbative l3{g) function up to two loops is 



2g 24 vr 



(llA^ - 2nf) + -— IsAN^ - lONuf - 3n ^ ^ 



4 TT V" ' "'^-^ A^ 



(36) 



with as = g {fi)/4:n. Eq.(^) is related to the vacuum energy through 



^^'"^) = \{Q,,). (37) 
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When we compare this last expression for (Q) with the one of Section 2 we can 
fix the behavior of the dynamical gluon mass as well as we can obtain a relation 
between the quark and gluon condensates. 

The calculation of the vacuum energy using the ansatze of the previous section 
is straightforward. However, it is more convenient to present the results of the 
calculation in different parts. The gluon contribution to the vacuum energy is 
obtained substituting Eq.p3|) into Eq. (]T5|) , and the result is 






where Hg is given by Eq.(|3J). 

Due to the form of the ansatz for the quark self-energy its contribution to the 
vacuum energy has to be computed only numerically. In the chiral limit we set 
rriQ = ruf = into Eq.(|T7p and easily verify that 

where h is given by 

h = z-^[ /. \'.^^! . 40 

V ln(r/VA^) ; 

When we consider bare massive fermions there is an essential modification in the 

calculation due to the mass renormalization. We initially consider the case where 

m/ << 1] and work with the ansatz of Eq. (^Oj) . Expanding the high energy part in 

powers of 1/P^, we verify that Eq.(p!7D reduces to 



(Qf) 



N 



^f 87r2 



dz z In 



z + ij] + nif) ri{ri + rrif) 



\ z + m'j z + 1]'^ 






7m-l 



hrf /z 



+ /^^ ...,„(! + VA)-^^ 
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8^ 



dzz (in^^^^^^^''^ ^^^ + ^fy 



,2 ^ _L ^2 



\ z + mj z + 7]^ 



-\ ■ '"" 



I dzz (ln(l + hr]^/z) 



(41) 



7m V V l + hrf/z^ 

where we kept only the leading term in mj in the ultraviolet part as well as the 
terms independent of the bare mass. A numerical evaluation of the above expression 
shows that (fi/);j- has a very small variation with nif up to masses of the order of 
100 MeV when we assume rj ~ 230 MeV with the minimum becoming 5% deeper 
for masses of this order. 

For very heavy fermions (rri/ » rj) the Schwinger-Dyson equations tell us 
that jTsI , |19] S(P^) K. mf{P'^), and it is not difficult to see that the fermion mass 
effect is basically erased from the vacuum energy Eq. (|l^) apart from a term coming 
from the mass renormalization, which survives no matter how small is the contri- 
bution from the dynamical mass. Therefore, for each heavy fermion of mass rrif we 
obtain 



N (—"(pip) mfiv) 

<"'>'./ --^i^^^t — ■ <''' 

Note that the renormalization group invariant quantity mj l—ijjipj appears after the 
mass renormalization. For light fermions the terms proportional to the bare mass in 
(ytf)if are significantly smaller in comparison with the terms due to the dynamical 
mass. For heavy fermions, i.e. for masses larger than a few GeV , all terms in 
{^f)f^f other than nif (—ipip) can be neglected without affecting considerably our 
forthcoming numerical estimates. 

To determine the parameter x which separates the perturbative and nonpertur- 
bative regions of the gluon self-energy and to find a relation between the condensates, 
we consider the equality 

(n(*'-)) = (fi,) + (fi;). (43) 

We fix the parameter x by working with the theory without fermions. We obtain 
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the following expression: 

8^ ^ ^ / 32^2 ^9^ 2 3 ^^ ' 

^ + X^)lnx+(^^)ln(x^ + l)). (44) 

Substituting rrig given by Eq.(0) and considering the P{g) function up to one loop 
we obtain 

(^-^ln(x + l)-Q + X^)lnx+(^^)ln(x^ + l))=^, (45) 

which gives 

X ^ 0.966797. (46) 



This value, being so close to 1, means that the gluon mass given by Eq.(|5^ is 
a reasonable scale to separate the perturbative and nonperturbative behaviors of 
the gluon propagator. In the following we assume a complete cancelation of (Qg) 
with the part of (i7*^*'"''\ coming from the gluonic contribution to the P{g) function 
when we use the above value of x- The phenomenological value of the gluon mass 
obtained in this way is also consistent with previous estimates 0. It is important 
to recall that in the OPE calculation of the gluon polarization tensor fermions do 
not contribute at leading order to the term proportional to (^G^^G^u\ and this 
cancelation is indeed complete. 

We can now consider the presence of fermions. The case of heavy quarks is the 
simplest one. Using Eqs.(|37|), (pSD, and (|42|) and taking into account the cancelation 
mentioned above, for Uf heavy fermions of equal mass we obtain 



i ^^G-G.„; - ^ (-^^) ,1 . -,„.) ^ "^-^ , (47) 



^f r^s^avr^ \ ^f^f 



im 



entailing 



rrif {-tfjtlj) ^ - l^^G^-G^,^ . (4J 
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Neglecting the second term of the [3 function, what is totally acceptable because a^ 
is small on the scale of heavy quarks, we obtain 

'^=12fl+7n^-7^4^V (49) 

which gives k k. 11.9, 12.6, 13.2, respectively for one, two and three heavy quark 
flavors. Thus, by calculating the vacuum energy as a function of the quark and 
gluon condensates we have reproduced the old result of Ref. fiy]. Note that this 
result appears naturally due to the cancelation between the gluonic parts of the 
vacuum energy and the fact that we can neglect the heavy quark contribution to 
the gluonic part of the effective potential, because the heavy quark contribution to 
the gluon mass proportional to the gluon condensate appears in the next order of as- 
We would like to note that the nonleading terms neglected in the operator product 
expansion, for example, the gluon condensate contribution to the quark self-energy 
S oc rrif (^G^'^G^yj jP^-, and vice versa the quark condensate contribution to the 
gluon polarization 11 oc mf \—'\\)-\\y) / P"^ may slightly modify our results, but even 
if we were intended to introduce them in the calculation we should remember that 
we do not have a reliable phenomenology to correctly describe their behavior at low 
momenta. 

The identity between the two calculations of the vacuum energy can also provide 
some information on the gluon propagator in the case of massless fermions. The 
physically interesting case is the one where the bare mass is zero but the chiral 
symmetry is broken dynamically. It is known that chiral symmetry breaking happens 
when the gauge coupling exceeds a critical value 



In principle this value (or a larger one |2^, ^) should be introduced in the I3{g) 



in Eq.(^). Unfortunately we do not know the behavior of the (3 function in this 
region! In this case we will use Eq.(^) at the extreme of its domain of validity, 
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which is the best that we can do at the moment. However, we still have another 
problem because our ansatz for the gluon polarization is valid only at leading order 
in Us (we do not explicitly take into account the fermionic contribution to the gluon 
polarization tensor, while we need to compare it to the two-loop (3 function). We 
proceed in a different way to the case of heavy fermions. For heavy fermions we 
used the identity between two ways of calculating the vacuum energy to obtain a 
relation between the condensates. Here, vice versa, we use it to obtain information 
on the ansatz for the gluon polarization. In what follows we modify our ansatz for 
the gluon mass performing the following replacement in Eg. ([3^) : 



ni 



9 a 



m Jl + agnfc), (51) 



where c is another parameter to be determined through Eg . (^31) . The parametriza- 
tion in Eg.(|51|) was chosen because the fermion contribution to the gluon mass 
proportional to the gluon condensate should depend on nf and as 

In the chiral limit we have to compute numerically the integral in Eg . (|39|) . As- 
suming that rj and Ag may vary around the values 0.2 and 0.3 GeV, we verified that 
the integral in this eguation contributes at most 10% of the total value between 
the brakets of Eg . (|39D , i.e. only the low energy region contributes effectively to the 
vacuum energy. If we neglect the high energy part of (Qf) and take into account the 
cancellation already discussed for the gluonic contributions in Eg.(^) we obtain for 
the chiral broken phase 



96 \7r ^7 V ' 47r V ' ^ N 



327r2 V 2 3"""'^' ■ '' V6 
-^lln(x +l)l/i,(l + a.n^c) ^^ 



+ I -^^^ I 1- (X^ + 1) 1 /^:(1 + "^%^) - .2^ ■ (52) 



With nj = 1, a, ~ a„ x given by Eg.(|6D, l^C^G^^ ~ (0.01) GeV^ fig, and 



-ipipj ~ (0.23 GeVY 1^, we obtain c ~ 8.7. This large value of c means that 
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the next order term in Eq.(|3^) in the case of hght fermions cannot be neglected. It 
would be interesting to compute the next order contribution to the dynamical gluon 
mass through the OPE including the effect of fermions and see how it compares to 
our estimate. Another numerical relation between the quark and gluon condensates 
could also be obtained for the case where the quarks have a small bare mass, using 
the values of x ^i-nd c determined above and Eq. (^11) for the fermionic contribution 
to the effective potential. However, such relation will be useful only with a bet- 
ter knowledge of the nonperturbative (3 function and of the next-to-leading order 
determination of the gluon mass through the OPE. 

5 Conclusions 

Starting from the effective action for composite operators [Q we determined the 
QCD vacuum energy as a function of the dynamical quark and gluon masses. The 
ultraviolet behavior of these masses can be related to the quark and gluon conden- 
sates through the operator product expansion, and the same is expected for their 
infrared parts, although, concerning this low energy behavior, we do have to rely on 
phenomenological data. With this procedure we obtain the vacuum energy described 
as a function of the condensates and bare quark masses. Equalizing this expression 
to the one obtained from the trace of the energy momentum tensor we were able 
to fix the behavior of the dynamical gluon mass. We verified that the gluon polar- 
ization tensor is fairly well represented by the ansatz we have chosen (Eq. (|33|) ) and 
that the gluon mass given by Eq. (|3^) is a reasonable scale to separate perturbative 
and nonperturbative behaviors of the gluon propagator. It is remarkable that our 
simple ansatz for the dynamical gluon mass gives such a good result. 

For heavy fermions we reproduced the old result of QCD sum rules mj (—tptp) ~ 
Yj-g (^G^'^Gi^iu) |]10[ which in our case comes naturally from the effective potential 



calculation. Note that this formula has been derived at one-loop and is essentially 



similar to the result of Shifman et al. in the sense that it is a consequence of the 
behavior of two-point functions. However, the concept of dynamical gluon mass is 
relevant here, without it we cannot compute an effective potential free of ambiguities 
(as a spurious infrared cutoff). We also discussed how the ansatz for the dynamical 
gluon mass could be improved and how it can be tested with the next-to-leading 
order calculation of the OPE. 

The method proposed here can be improved in several ways. We certainly need 
a description of the /3 function in the nonperturbative region, without this we do not 
believe that we can go far away with these calculations in the light fermion domain. 
The other improvements are related to the approximations that we have performed 
to compute (Q), i.e. the choice of the ansatze for quark and gluon propagators. The 
computation of (Qg) has shown to be quite reliable |]14|, and in spite of neglecting 



the contribution of diagrams with quadrilinear vertices and ghosts, the final result 
is similar to the Hartree approximation made by Cornwall in his calculation of the 
vacuum energy [0, and both calculations reduce the two-loop contributions to a 
term whose main contribution to the effective potential is proportional to 11^. In 
a more rigorous calculation the two-loop effective potential should be calculated 
including all diagrams, as well as the full vertices because they are fundamental 
for cancelations present in the full calculation |2^. As we computed the effective 
potential at the minimum of energy we escaped from this very difficult task, and 
also minimized the effect of the approximations p!T |. 



Another point that could be discussed at length is the question of the dynamical 
masses ansatze. It is obvious that the desirable scenario would be a full numerical 
solution of the Schwinger-Dyson equations concomitantly with the effective poten- 
tial, but this is also an enormous work. The ansatze we used are in agreement 
with the phenomenology at low momenta and the OPE at high momenta. Notice 
that the important point is the connection between the dynamical masses and the 
condensates, and it is not relevant what mechanism generates the condensates or 
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masses as long as the relation predicted by the OPE between the condensates and 
the masses holds. 

Finally, the idea of the calculation is very simple, and we were able to recover, 
within the approximations of this method, the relation between quark and gluon 



condensates for heavy quarks obtained by Shifman et al. many years ago |^, and 
verify that the vacuum energy is well described by a dynamically massive gluon, 
whose momentum dependence is displayed in Eq.(^) and Eq. (|3^) . 
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Figure Captions 

Fig. 1 Two-particle irreducible vacuum diagrams contributing to the effective 
potential. 

Fig. 2 Diagrams contributing to the gluon polarization tensor. 
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